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We study various representations of infrared effective theory of SU{2) glu- 
C<") , odynamics starting from the Abelian monopole action derived recently by 
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numerical calculations in the Maximal Abelian projection. In particular we 



o 

(N 

derive the string model and the dual Abelian-Higgs (dual Ginzburg-Landau) 



model which corresponds to infrared SU{2) gluodynamics. It occurs that the 
classical string tension in the string model is of the same order as the string 



<—| ! tension in quantum SU(2) gluodynamics. 
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I. INTRODUCTION 

The knowledge of an infrared effective theory of QCD would be very important for the 
analytical description of hadron physics. Before the derivation of a low-energy effective 
theory of QCD we have to explain the most important non-perturbative phenomenon, the 
quark confinement. Wilson's lattice formulation |lj shows that the confinement is a property 
of a non- Abelian gauge theory of strong interactions. At strong coupling the confinement is 
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proven analytically. At weak coupling (corresponding to the continuum limit) there are a lot 
of numerical calculations showing the confinement of color. The mechanism of confinement 
is still not well understood. One of approaches to the confinement problem is based on 
the searching of the relevant dynamical variables and construction of an effective theory for 
these variables. 

From this point of view the idea proposed by 't Hooft || is very promising. The pro- 
posal is based on the fact that after a partial gauge fixing (Abelian projection) the SU(N) 
gauge theory is reduced to an Abelian [/(l)^ -1 theory with N — 1 different types of abelian 
monopoles. Then the confinement of quarks can be explained as the dual Meissner effect 
which is due to the condensation of these monopoles. The QCD vacuum is dual to the ordi- 
nary superconductor: the monopoles playing the role of the Cooper pairs. The confinement 
occurs due to the formation of a string with the electric flux between quark and anti-quark, 
this string is a dual analogue of the Abrikosov string [[|. The described mechanism of 
confinement is usually called the dual superconductor mechanism. 

There are many ways to perform the Abelian projection, but only in the Maximal Abelian 
(MA) gauge ||] many numerical results support the dual superconductor picture of confine- 
ment |J (see, for example, reviews @|7|). These results suggest that the Abelian monopoles 
which appear after the Abelian projection of QCD, are relevant dynamical degrees of freedom 
in the infrared region. 

In compact QED, due to the periodicity of the gauge field, there exist monopole excita- 
tions which are responsible for "charge confinement" ||. The effective theory in terms of 
monopoles was studied both analytically || and numerically [[K]]. It was found [10| that 



the monopole action contains all information about confinement in this theory. The role of 
monopoles in the confinement in non-Abelian gauge theories is also non-trivial. The well 



known example is confinement in the 3-dimensional Georgi-Glashow model J] 1]] which is due 
to the 't Hooft-Polyakov monopole fl"2|| . 

The effective monopole action for the MA projection of SU(2) gluodynamics was obtained 
by Shiba and Suzuki fH| using the generalized Swendsen method [[H|]. Assuming that the 



action contains only terms which are quadratic in monopole currents, they found that the 
action has the form theoretically predicted by Smit and van der Sijs [|14]]. However the 
considered physical distances are rather small and the main assumption, the quadratic form 
of the monopole interactions, is not well justified. 

The monopole action which contains only quadratic terms corresponds to the dual 
Abelian Higgs theory in the London limit, the lagrangian of which is 

L dual (A, $) = -L(d [/X A,]) 2 + \\{d, + iA,)<$>\ 2 + A(|$| 2 - V 2 ) 2 , (1) 

where A^ is the dual gauge field and $ is the monopole field. In the London limit, A — ► +00, 
the radial part of the monopole field is frozen, |<&| = 77. As we discuss in Section || the 



corresponding monopole action in the lattice regularization is ||T4| . |T5 

S mon (*k) = Y(* k > A ~ 1 * k ) + P\\* k W 2 ' ( 2 ) 

where the first term in the r.h.s. corresponds to the Coulomb-interaction of monopole 
currents *k, A -1 is the inverse lattice Laplacian. The last term in (Q) is proportional to the 
length of the monopole trajectory. In (|2|) and below we use the notations of the formalism 
of the differential forms on the lattice |16|, for the brief explanation of the notations see 



Appendix A of the review 0. 

It is not necessary that SU{2) gluodynamics in the abelian projection corresponds to 
the dual Abelian Higgs theory (JJ) in the London limit (A — > +00). The finite value of A 
was firstly discussed in ref. |17j where the dual Ginzburg-Landau (DGL) theory was studied 



as an effective theory of QCD. From the analysis of the experimental data it was found 
that DGL is near the border of type II and type I superconductor (Bogomol'ny limit). The 
similar conclusion has been made in Ref. |LJ| where a non-abelian dual theory of QCD has 



been considered. At present there are a lot of numerical data JT9[ showing that the lattice 
gluodynamics is in a sense equivalent at large distances to the dual Abelian Higgs model 
near the Bogomol'ny limit. In the present publication we present new facts in favor of this 
statement. 



Below we show that for finite coupling A the additional terms in the monopole action (0) 
appear. This is consistent with the recent numerical analysis of the monopole interactions 
on the large (48 4 ) lattice performed in ref. J20 |; where it was shown that the effective 



monopole action contains not only quadratic interaction but also 4- and 6-point couplings 
of the monopole currents. The action seems to satisfy the scaling behavior, that is, it 
depends on the physical length b = na{j3) alone, where n is the number of the blocking 
transformations and a(/3) is the lattice spacing. This scaling behavior is consistent with the 
behavior of the perfect action on the renormalized trajectory. In Section |TTT| , we analyze the 
new numerical results on the monopole action. The couplings of the monopole action are 
determined by Swendsen's method from the monopole current configurations extracted in 
the MA projection from SU{2) lattice gauge fields. 

In Section [IV] we derive the string action from the monopole action. We thus obtain 
the effective string action for lattice gluodynamics. There exist many attempts to relate the 
string theory with a local field theory. Nielsen and Olesen |2T] have found quantized magnetic 
flux lines as static solutions in the Abelian Higgs model, they also mentioned that strings 
behave similarly to the Nambu string WM. Forster et al. [E3] have proved that, in the zero- 



width limit, flux lines move like Nambu strings. Using the collective coordinate formalism 
p3^p5| they have constructed an effective string model. There are many speculations about 
the possible relation between string theory and non- Abelian field theory [f26|j . 

In the London limit (A — > oo) the Abelian Higgs model, can be exactly rewritten as 
the string theory with a non-local action [[25|j27| -|29H . Thus in the London limit the infrared 



effective theory of QCD developed in Ref. [17| can be explicitly rewritten as a string theory 



in the London limit. This fact is natural since the dual picture of QCD vacuum suggests the 
formation of QCD strings. The formation of string between quark and anti-quark leads to 
the confinement phenomenon. The effective action for the QCD string contains in addition 



to the Nambu-Goto term the so-called rigidity term p7|j30|j3T[] and other more complicated 
interactions. 

The main task of this publication is to interrelate various models (monopole current 



model, string model, Abelian Higgs model and others) which are equivalent to lattice glu- 
odynamics. In Section [II] we briefly discuss various models which are equivalent to the 
quadratic monopole action (|2]). In Section [Tfl] we calculate the couplings for the monopole 
action which corresponds to SU(2) gluodynamics in the MA projection. In Section [TV] we 
discuss the effective string action corresponding to lattice SU{2) gluodynamics. In Section [V] 
we discuss the DGL theory (an effective dual Abelian-Higgs model) which is equivalent to 



the theory with the monopole action described in Section |Tj. We study this effective model 



in three different regions of coupling parameters of the DGL theory. The details of the 



calculations are collected in Appendices A-B. The conclusions are given in Section VI. 



II. MONOPOLE ACTION IN THE LONDON LIMIT 

In this Section we discuss several representations for the theory with the quadratic 
monopole action (0). The basic method we use is a generalization of the transformation 
of Berezinski |[32|| , Kosterlitz and Thouless |33[ and the duality transformation originally 



suggested by Kramers and Wannier [53] for the 2-dimensional Ising model. The Berezinski- 
Kosterlitz-Thouless (BKT) transformation relates the two-dimensional XY model and the 



"vortex + spinwave" system |32|.|33|.|3l| . Banks et. al extended this method to the system 
with Abelian gauge symmetry, and found a "monopole + photon" representation of lat- 
tice gauge theory in 3-dimensions and a "monopole current + photon" representation in 
4-dimensions @. The monopole partition function corresponding to monopole action (0) 
can be transformed into the partition function of the dual Abelian Higgs model @,|l^] and 



into the partition function for the string theory of hadrons on the lattice |2S 
The partition function in the case of quadratic monopole action is: 

z mon = J2 eX P (S mon (*k)) , (3) 

*fc(* ci )g^ 
S*k=0 



where the monopole action S mon is defined byQ eq.(H). The condition S*k = means that we 
sum over all closed (conserved) currents on the dual lattice. Applying the inverse duality 
transformation (see | Appendix A| ) to eq.(|3|) we get the partition function: 

z au = f D * A {* Cl ) f D* V {%) T exp{--^-||rfM|| 2 -^||rfV + 2vr1- Mf) . (4) 

Now we show that (|j) is the lattice version of the (dual) Abelian Higgs model in the 
London limit. The direct lattice analogue of the action of the Abelian Higgs model ([I]) is 



+A E ($*(x)$(x) - l) 2 + E $*(x)${x), (5) 

X X 

where U = exp (iA^) is the dual gauge field, dA = A P is the field strength tensor and 
$(x) = p^exp (iip x ) is the complex Higgs field. 
Action (|5]) can be rewritten as follows 

S A - u -[A,p,ip] = -L \\dA\\ 2 - 2^ p xPx+ii cos {(dtp) Xtfi + A,} 

Urn x,(i 

+AE^-i) 2 + E^- (6) 



It is convenient to modify the action using the Villain formulation by replacing 
exp {a cos ip} — > J2i<=% ex P { a ~ f ("0 + 2vr/) k and the partition function of the Abelian 
Higgs model becomes: 

+ 00 +7T +00 



-|-uu 1-71 T"U • 

Z A - R -= I VAJV V /p/^exp -- T ||dA|| 2 - 7 EP^ + A(^ + (^ + 27r/ 

^ L 9m x ^ 

7E(P^-P^A) 2 -E((l-87)P 2 +A(p^-1) 2 )}. (7) 

X,[l X ) 



x In this Section we restrict ourselves to the form of the quadratic action given by eq.(^), however 
the same transformations can be performed for any quadratic monopole action. 



The continuum limit of the Villain action is still the action defined by eq. ([!]). In the London 
limit A — > oo, p x = 1 and the partition function (|7|) is equivalent to the partition function 
(||). The relation of the lattice lagrangian which enters the partition function ([|) and the 
continuum lagrangian (|1|) of the abelian Higgs model is obvious. d*A is an analog of the dual 
field strength tensor dht s A u ] and the term \\d *<p + 2n *l — *A\\ 2 corresponds to | (d^ + iA^Q>\ . 
The summation over I appears due to the compactness of the phase ip of the Higgs field 
$ = |$| e l{p . The radial component |<&| of the Higgs field is frozen, since A — > +00. 

Performing the BKT transformation for the field *ip(*Co) in eq.(|J) we get the string 



representation of the model (see | Appendix A| for details) [29] 



7T 2 / / . . o\-l 



o-(c 2 )G^,(5cr=0 k P 



Z StT = E exp -- a, A + Mi a , (8) 



where 



9v 



M l = %> (9) 

is the mass of the dual gauge field A. The integer-valued 2-form a represents the closed 
world surface of the string formed by the electric flux. The propagator (A + M\) in eq. ([8]) 
means that the string elements interact via the exchange of the massive gauge boson. The 
condition 5a = means that the world sheets of Nielsen-Olesen strings must form closed 
surfaces. This interpretation of a is clear from the following considerations. Consider in 
representation (0J) the expectation value of the 't Hooft loop operator which introduces the 
quark-ant i- quark pair: 

OO 7T 

< H ( c ) > = zam. J D * A ( * c / D V( %) E (io) 



-00 



H(* Cl )<=% 



exp { \r\\d*A + 2n*S c \\ 2 - — lid V + 2tt 1 - Mil 2 

1 ^ f 7T 2 / /. .,,\-l \ 27T 2 



£ exp _^( <7 ,(a + JW1)- 1 (t)-^(c,(A + ^)- 1 c) , (11) 



tr(c 2 )ez: (. f K ' ™ 

Sa=C 



here C is the quark-antiquark trajectory; Sq is an arbitrary surface spanned on the loop C: 
SSc = C. The summation in ([11]) is over all closed surfaces (virtual glueballs). Effectively 



it means that we sum also over all surfaces spaned on the contour C. Thus a represents the 
worldsheet of the color electric flux. 

Finally we mention that the abelian monopoles in model (|3]) can be considered as defects 
in a compact gauge theory similarly to the compact QED. The corresponding "compact 
QED" representation of the partition function (^) is: 

+ 7T 2 

Z compact = / DO Y exp{-^||d# + 27m|| 2 --^-(d# + 27rn,A(d# + 27rn))}, (12) 
J txr-v, L on 2 An 2 > 

where n is the integer- valued two-form and the action is invariant under U(l) gauge trans- 
formations, 9 — ► (9 + da) mod27r. Applying the standard BKT transformation [[],|3^,[£| to 
the field 9,9 + lixn = 9 nc _ + 2irA~ 1 5k, (*k is the closed monopole trajectory, 5*k = 0, and 
9 n _ c _ is a non-compact field), we get the partition function (|j) from z compa,ct (|T2| ). 

III. THE MONOPOLE ACTION FROM GLUODYNAMICS 

In this Section we present the results of the numerical investigation of the effective 
monopole action in SU(2) gluodynamics. We generate the SU{2) gauge fields on 48 4 lattice, 
perform the Maximal Abelian gauge fixing and extract monopole currents. 

As we discuss in the Introduction, there are strong indications |TT|-|T9|j, that the gauge 



boson mass is of the order of the Higgs boson mass in the effective Abelian Higgs model 
which corresponds to the infrared lattice SU(2) gluodynamics. Thus A is finite in the 
effective action ([I]) and (](]). As we show in Section |V| for finite A the monopole action 
contains iV-point (N = 2, 4, 6, . . .) interactions. Hence we have fitted the numerical data for 
the monopole currents by the partition function^]: 

Z mon = Yl exp{-S mon [*£;]} , (13) 

S*k=0 

where the action <S mon [*fc] contains 2-,4- and 6- point interactions of the monopole currents: 



'For the details of numerical calculations see Ref. [2C] 



S™»[*k] = ^ (% A~ l *k) +p(b) \\*k\\ 2 -g(6)E f £ **£/. 
+ r(6)E E *CC( 3, ^ + *0+E/.(^?[1. (14) 

£ /x,i/,<5=— 4 i 

The monopole currents *fc are defined on the dual lattice. They are integer valued and 
conserved (<5*fc = 0). S\ [k] are additional 2-point interactions [IIJ of monopole currents 
which are introduced to check whether there are any corrections to the Coulomb interaction 
(*k, A _1 *A;). Our numerical results show that these corrections are negligibly small in the 
infrared region, and we neglect them in what follows. 

To determine coefficients g^b), p(b), q(b), r{b) numerically we use the modified Swend- 



sen method [10,T3[ (which is a kind of inverse Monte Carlo method) for monopole currents 
in the MA projection of lattice SU(2) gluodynamics [Q. In order to study the renor- 
malization group evolution of the effective action in the infrared region we consider the 
n-fold extended type-II monopole currents |3jJ (n = 1,2,3,4,6,8) which emerge from the 
block-spin transformation on the dual lattice. 

It is interesting that the couplings of the monopole action ( |14"D seem to depend only on 
the physical length b = n- a(/3) when n is larger than 3. This fact means that we are working 
on the renormarized trajectory. In order to express b in terms of the physical (dimensional) 
units we measure the dimensionless string tension k{(3) (in lattice units) pi| . Then a (j3) = 



K( / 9)//€phy S , where K p h ys is the physical string tension (which is about 440 MeV in SU(3) 
QCD). 

The best fit for the renormalized couplings to numerical data for g^ib), p(b), q{b), r(b) 



gives 3 : 



9l(b) norwo ^/4vr\ 2 



).300(22) ( 



\9(b). 



3 We use numerical data for g^b), p(b), q(b) and r(b) extrapolated to n — > oo limit for each value 
of b. This extrapolation corresponds to continuum limit of lattice theory, see Ref. [20] for details. 
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where 



Pip) 
q(b) 
r{b) 



9 (b) - 

A = 0.88(7). 



11 , L i 

In 1 + 



(bAY 



0.527(25) 

£3.53(7) ' 

0.0178(15) 

66-18(13) ' 

0.83(8) x 10~ 4 
510.05(17) • 



17 , 

447T 2 



In 1 + 



(bA)' 



(15) 



(16) 
(17) 



These fits are shown in Figures [T|-^]. The 6-dependence(16) of g is consistent with the 2-loop 
running coupling for the small b region and reproduces the experimental (numerical) power 
behavior for the large b region. 



IV. THE LATTICE QCD STRING 



-1/2 



It follows from Figure || that the 6-point interaction is negligibly small for b > 1.0 
We neglect it for simplicity in what follows, since we are interesting in the infrared effective 
theory which corresponds to large b values. The introduction of the small 6-point interaction 
is not difficult. Our starting point is the monopole partition function and we discuss the 
string representation for the monopole action (cf. eq.fTTJ])): 

£ e W \- 9 ^(*k,A-^k)- P \rk\\ 2 + q ^(j:*kl\}. (18) 



Z 1 



"fc(*ci)e^ 

S*k-0 



K ix=-A 



The partition function ( |18"D is divergent at large densities of the monopole currents *k 
due to positivity of the coefficient q. However, the partition function is a part of the total 
partition function with n-point, n > 6, being omitted. Due to the higher-point interactions 
the total partition function is convergent and the expectation value of the monopole density 
is finite. Therefore at sufficiently small values of q we can treat the coupling q in eq . (|18D 
perturbatively: 



10 



*fe(*ci)e^ xx \a=-4 / 

<S*fc=0 

= II(l+«(Ej!-) +Otf))2[B]| M , (19) 

where the generating partition function is: 

Z[B]= E exp{-^(*k,A- l *k)-p\\*k\\ 2 + (ek,*B)}. (20) 

<5*fc=0 



Using the transformations of |Appendix A| we get: 



Z[£] = V expj-— (Wa + cLB, (A + M|) * (2trx + cLB 
Substituting this expression in eq.(|T9"D we get the following partition function 

Z mon ocZ stT = J2 exp{-S (0 \cx)-qS {1 \cx) + O{q 2 )}, (21) 

o-(c 2 )£^,5o-=0 

where 



7T 2 , / . . , \-l 



5 (0) (a) = ^-(a,(A + Ml)" a), (22) 

S (1) (tr) = E E (2T«,^^W. + 4^,^,^^^ - 4^, J , (23) 

x n,v=— 4 



^ = - (A + Mjp^ (24) 

py 2p v J x,v;v,v 

By a straightforward calculation one can show that 

T^ ;x ,, = |-(ZW0, 0, 0, 0) - P A f A (l, 0, 0, 0)) ( for /x = i/ ) 

= ^(2P Ma (1,0,0,0)-Pm a (0,0,0,0)-P Ma (1, 1,0,0)) ( for // ^ z/ ) , (25) 
= ( for /i = — z/ ) , 

where Ma is effective dual gauge boson mass (^) and T> Ma [x) = (A + Mj|) (x) is the 
massive propagator. Combining eqs.(PI|)-(|2*5D we get the string action: 
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S s >] = -(a,V MA a) + ^1(V Ma (0, 0,0,0) 
p p A 



-V MA (l,0,0,0))-Y,J2 ( V mJv) 2 +?*L(2D Ma (1, 0,0,0) (26) 

x M =-4 ^ P 

-v Ma (o,o,o,o)-v Ma (i,i,o,o))-Y: E K^.i^,^),, 

E(E>«>C) 2 + otf 



(MtM 



The leading part of the action (the first term in the r.h.s. of fl26|)) comes from the 
self- interaction and from the Coulomb interaction of the monopole currents (see Section |J); 
there are also corrections due to *fc 4 -interaction. Action ( pB"D is calculated as an expansion 
over the parameter q which is small in the large b region. It is interesting to note that the 
4-point interaction in the monopole action induces the 4-interaction in the string action (last 



term in eq.fl26|)). The presence of the 4-point term is an indication of the finiteness of the 
dual Higgs boson. Thus the local four-point monopole interaction leads to the non-local 
four-point string interaction. 

The classical string tension is defined as follows: 

S(a flat ) = K th ■ Area(a flat ) , (27) 

where <r flat is an infinitely large flat surface. 
There are three contributions to n t h_'- 

^ = ^ ) +g(^ ) ' // + 4 ) ' /y ) + 0(g 2 ), (28) 

where n t ^ is the contribution from the leading term S^°\ eq.(^), and n th , n ti ^ are 
the contributions from, respectively, the two- and four-point interactions of the action S^\ 



Calculation of K th gives: 



K 



(0) _ ^TlP) 



-VWM, (29) 



th — _, ^M A 



P 

where T>m (x) = (A + M\) is the massive propagator in two space-time dimensions. 



M A l 
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The calculation of k[^ is more involved. Consider an infinitely large flat surface in the 
plane (3, 4) which is defined by equations X3 = £4 = 0. We note that for the flat surface 
*(J a p(x) takes non-zero value only when sitting on the (1, 2)-plane : 

V? 2 at (x 1; x 2 ,0,0) = -V 2 fl f(x 1 ,x 2 ,0,0) = 1 

V?f(x) = (for *, jV 1,2) (30) 

Then the two-point contribution to the string tension is: 

^) = 1 2 ^! ( P Ma(0)0)0)0) _P Ma(1)00)0)) .(P« (0,0) -2^(1,0)) 

8tt 2 



-(2V Ma (1 } 0, 0, 0) - V Ma (0, 0, 0, 0) - %(1, 1, 0, 0)) 



p 3 

(2V ( S A (1,0)-Vfl(0,0)-V ( S A (1,1)). (31) 



The four-point contribution to the string tension can be obtained from eqs. (|26| , |24|) : 

2 



t ( i ) '^ = -^E( E (^(x + n(«))-^(f+n(a)+ea) 2 ) , (32) 



r 2 

P x a/3=-2 

where n(a) = if a > 0, and ra(a) = —(1,0), if a = — 1 and n(a) = —(0, 1), if a = —2. 
Here (ea)/3 = £ a p, and e^ is the modified antisymmetric tensor in two dimensions: s a p = 
sign(a) sign(/3) e a/ g, en = e 2 2 = 0, 612 = — e 2 i = 1. For large fe-region the term /% is much 
smaller than n th and n th . 



For b > 1.0 



K 



-1/2 



the theoretical string tension is almost constant and reproduces 
the physical string tension unexpectedly wellQ, Kth ~ 1.3 (see Figure [5]). Note that K t h in 
eqs. (^,^,[T!],[^) is the string tension of the classical string, since it is simply the coefficient of 
the Nambu-Goto term in the string action (p6|). On the other hand it is close to the string 
tension of the quantum lattice SU(2) gluodynamics. This fact means that the quantum 
corrections to the string tension are small. 



4 In our normalization, if K t h = 1 then in physical units it equals exactly to the string tension of 
SU(2) gluodynamics obtained numerically. 
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V. RELATION OF THE MONOPOLE MODEL AND DUAL ABELIAN HIGGS 

MODEL 

From numerical calculations we know couplings (15) of the monopole action flT4|) . It 
would be interesting to derive the Abelian Higgs model which corresponds to flT4|) , (|15|) . In 
this Section we show how to solve the equivalent problem, i.e. how to derive the partition 
function of the monopole current from partition function of the lattice Abelian Higgs model 

(0)- 

Inserting the identity 



const. = Vet 2 (A^p x p 



X+fl) 



x 7 PFexp { - £ — ^ (^ - 2*7PxPx + a {A + {dip) + 2irl)f ) (33) 

into the partition function ([?[) and integrating over the fields A and tp, we get the following 
monopole representation of the partition function (cf. eqs.(3,P^)): 

2 oc 2 mon = J2 exp{-S mon [*k}}, (34) 

s*k=a 

where the summation is over all closed monopole trajectories *k. The monopole action 
S mon [*k] is: 

S mon [*k] = 5g lon [*k] + S£ on [*k] , (35) 

where the Coulomb part 

comes from the integration over the gauge field A. The Higgs part 

exp {-SJT [**]} = / ^P" 2 exp - 7 ||dpf "EE I ~ -V(p)\, (37) 

( x M=l 4 TPa;Pa;+At J 

includes the result of the integration over the modulus of the Higgs field |$| = p and V (p) 
is the Higgs potential: 
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^(p)=E{a(p^.-^) 2 + A(1-(^) 2 )}, (38) 

2 2A + 87 - 1 



2A 



The monopole action denned by eqs. (p5|) -(|38D is the main result of this Section. The 
integration over p in eq.(^) can not be performed exactly and in order to relate the monopole 
action (|35"D-(|3~8"D with the action (|I4] ) ,(|T5 D obtained from numerical calculations we have to 
develop the approximate methods of the evaluation of the integral ( |37|) . 

The simplest case is the London limit (A — > 00), the radial part of the Higgs field is fixed, 
p x — 1, integral ( ft?D is trivial and S mon in eq.(^) coincides with S mon given by eq.(0). This 
case corresponds to the effective action at very large distances, b — > 00 (see eqs.(|I~4^,(|T5|)). 
To get the effective infrared lagrangian of lattice gluodynamics at finite physical scale b we 
study the ranges of parameters A and 7 at which integral (^) can be calculated analytically. 
The regions are: I. A ^> 1, 7 ~ 1; II. A = cry 3> 1, c ~ 1; III. 7 <C 1, A ~ 1. The monopole 
actions for these regions are given in Appendix B. 

Our numerical calculations are restricted to the region b < 4, just in this region we get 
the couplings ([15]) by the fit to the numerical data. It occurs that regions I, II, III do not 
correspond to couplings ([15|) for b < 4. In other words, eqs. (p.2j ),( p.8| ) and ( |B.11|) have 
no real solutions for couplings A and 7 if we substitute p(b) and q(b) defined by eq. (|i~5|) . 
This result is natural, since we know from phenomenological fT7| , |T8"|| and quasiclassical (19 



analysis that for b < 4 SU(2) lattice gluodynamics corresponds to the Abelian Higgs model 
near Bogomol'ny limit 7 ~ A ~ 1. But in this region integral ( pTf ) can not be estimated 
analytically. We discuss regions I, II and III in Appendix B since in the future numerical 
calculations for large values of b the couplings A and 7 may lie in one of these regions. 

VI. CONCLUDING REMARKS 

In order to study the effective infrared action of lattice gluodynamics we performed the 
following steps. 
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1. The abelian monopole action is extracted from the SU(2) gauge fields in the Maximal 
Abelian projection. 

2. The couplings of the monopole action are calculated from the ensemble of the monopole 
currents using the modified Swendsen method pO|JT3|j20|] . It occurs that the couplings 



depend only on the physical length b = n- a(j3), thus we are working very close to the 
continuum limit. The coupling of the four-point interaction of the monopole currents 
*k 4 is definitely not zero, thus the corresponding dual Abelian-Higgs model is far from 
the London limit for the considered values of b (0.5fm < b < 2.5fm). 

3. We derive the relations between the couplings of the monopole action and the dual 
Abelian Higgs model near the London limit. These relations can be used to get the 
parameters of the effective dual Abelian Higgs model which will be obtained in the 
future numerical calculations of the monopole couplings on large lattices. 

4. From the effective monopole theory we derived the effective string theory for lattice 
gluodynamics. It occurs that the classical string tension of the effective string model is 
close to the quantum string tension of SU(2) lattice gluodynamics. Probably, it means 



that the (quasi-) classical string theory defined by action fl2"E| ) is a good approximation 
for infrared gluodynamics. 
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Appendix A. Transformations for Quadratic Abelian Higgs Action 

Using an auxiliary field *A( *ci), the partition function (|3|) can be written as 

Zoc 7 D*A(* Cl ) Yl exp{-yr\\d*A\\ 2 -i(*A, *k) - p\\ *k\\ 2 } ■ (A.l) 

Introducing the phase *<£>(* c o)> the current conservation law can be expressed by the delta 
function: 



5(5 *k) = / D*ipexp[i(*ip,5*k)]. (A.2) 

— 7T 

Substituting eq. (|A.2| ) in eq. (|A.l| ), we get 



OO TV 

Z<x J D*A(* Cl ) J DV(*co) £ 



— 7T 



*fc(*ei)e^ 



exp^ _||dM|r + i(dV- % ^)-p||^llf • ( A - 3 ) 



'in 



Our next task is to integrate out the monopole current *k. For this purpose, we replace 
the integer- valued field *k by a real- valued field. This manipulation is accomplished by the 
Poisson summation formula, 

[DFj2exp[2m(FJ)]f(F) = £ /(*), (A.4) 

where / is an arbitrary function. 

Applying this identity to eq.( |A.3|) we get 



OO OO 7T 

Zoc J D*A(* Cl ) J D*F(* Cl ) J £>V(*co) £ 



— oo — oo —TV H cijG-aJ 

expl ^-||rfM|| 2 + i(rfV + 27r1- M, *F) - pll *F|| 2 1 . (A.5) 

[ % J 

The Gaussian integral with respect to the auxiliary field *F(*ci) leads to eq.@. 
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To get the string representation from eq. (|4]) we should perform the BKT transformation 
for the field 1(*Ci). Substituting the decompositions 1 = *s[*a] +d*r, d*s = *a, Vn.c. = 
*ip c + 2 / k(5A~ 1 *a + V) to eq.(||) we get the partition function in the following form 



Zoc J D*A(\) J D*<p(*c ) £ 



-oo 



: a(*c 2 )eZ,d*<7=0 



exp(--^||rfM|| 2 -^||rfV + 2vr(5A- 1 V- *A\\ 2 ] (A.6) 

Choosing the gauge d*(p = and integrating over V(* c o) an d *A(*ci), we obtain lattice 
surface (string) model (|8]). 

Appendix B. Monopole Action in Three Regions 

In this Appendix we present the results of the exact calculation of S^ on for three different 
regions of the parameters 7 and A. The details of this calculation will be published elsewhere. 

Region I: A ^> 1, 7 ~ 1 
In this region we integrate in eq.([Tj]) over the radial mode p using the saddle point 
expansion method since the coupling A is large. 
The result is: 

ST n i*^=Pl\*k\\ 2 -qj:(j:*kl l ) +0(A- 2 ), (B.l) 

x W=4 / 



where 



' +0(A~ 2 ). 



2 8 7 2 A 



Region II: A = 07, c ~ 1 

To compute integral ( fi7\) in this region we can use the saddle-point method with respect 
to 7 and A. 



In the leading order the monopole action is: 

3 



^ = E^ } A^ + 0(A- 4 ), (B.3) 



where 



7 (i) c 



2 



s i" = m^T^ '; (B.4) 

^ = 8^11*11"= < R5 > 

S<? ) = --^-j(p,(A + 4(c + 4)r 1 p). (B.6) 

(c + 4) v 7 

(B.7) 

Taking into account only local contribution to the monopole action we get the action of 
form (|B.1|): 



c 2 



* = 4Mt+T) + °W (R8) 

^ 0i0 (4(c + 4)) + O(A- 4 ). 



2 6 A 3 (c + 4)^ 
where T> xy (m) is the propagator for scalar particle with the mass m: (A+m 2 )'D Xjy (m) = 5 x ^ y . 

Region III: A < 1, A ~ 1 

The saddle-point approach does not work in this region. And we use the perturbative 
expansion in powers of 7. The results is: 

s%r = ^{k,A- i k)+ P \\k\\ 2 -q'E( e *£j 

Z x M =-4 

+-E( E *£j + *E( E *£j +o( 7 ), (b.9) 

£ |U=-4 X f_L=~A 

where the coefficients are given by the following expressions: 

1 1 . / / -25 % , A o , 4u 3 m 4 

p = — In — hp, p= h 4 «! — 3 m 2 



7 12 3 4 

1 u + 104 ui - 114 u 2 + 56 m 3 - 1 

96 

-5 mq + 18 Mi — 24 w 2 + 14 M3 — 3 W4 



-35 u + 104 ui - 114 u 2 + 56u 3 - 11 u 4 

<Z = gg , (B.10) 



96 

M — 4 «! + 6 U 2 — 4 M 3 + M 4 

384 ' 



(B.ll) 
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the functions U{ are: 

MA) = -ln^(A), fc = 0,...,4, (B.12) 

and the functions g/i(A) are: 

^i(A) = ^(e- A + (-l + 2A)^ (A)), 

^ (A) = (2X? ((_1 + 2A) 6 ~ A + (1 " 2A + 4A2) ^° (A)) ' 

^( A ) = T^T^ 1 + 4A ') e " A + (- 1 + 8A3 ) ^( A )) ' ( B - 13 ) 

[ZA) 

04(A) = T^i((- 1 - 4A + 8A ' + 8A ') e ' X + (1 + 4A - 12A 2 + 16A 3 + 16A 4 ) g (X)) , 
(2A) 

where Erf(x) is the error function: 

+00 
Erf(x) = 4= / e- t2 dt. (B.14) 
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FIGURES 
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FIG. 1. The numerical data and fit ( |l5| ) for the coupling of the monopole action g^ib). 
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FIG. 2. The same as in Figure IT] but for p(b). 
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FIG. 3. The same as in Figure 111 but for q(b). 
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FIG. 4. The same as in Figure [T] but for r(6). 
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FIG. 5. The theoretical string tension eq.(28,2£,31,32) and numerical data for the full SU(2) 
string tension vs. b. 
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